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In this paper we introduce a novel method for fabricating MgB2 /MgO multilayers and demonstrate 
the potential for using them as a new model for disordered superconductors. In this approach we 
control the annealing of the MgB2 to yield an interesting new class of disordered (amorphous) 
superconductors with relatively high transition temperatures. The multilayers appear to exhibit 
quasi-two-dimensional superconductivity with controlled anisotropy. We discuss the properties of 
the multilayers as the thickness of the components of the bilayers vary. 

PACS numbers: 74.70. Ad, 74.25. Op, 74.62. Bf, 78.40.Pg 



I. INTRODUCTION 

Magnesium diboride (MgB2 ) has proven to be a dif- 
ficult material to synthesize in thin film formJ^ The pri- 
mary source of this difficulty lies in the very high vapor 
pressure of Mg under the conditions of growth. Similar 
to single crystals, the best method to fabricate high qual- 
ity thin films of MgB2 has proven to be the Mg diffusion 
method, which yields results that are almost independent 
of the substrate.^ A general approach to the fabrication 
consists of the deposition of a disordered precursor film, 
usually with excess Mg, followed by an anneal to form 
crystalline MgB2. This method will yield high quality 
films with high transition temperatures. Such films may 
be important for future device applications such as tun- 
neling devices that work in the range of 25 K, and may 
have the potential to replace low transition-temperature 
materials (NbTi, NbaSn) in high- field applications. 

Besides its very high transition temperature (~ 40 K), 
its low anisotropy (< 6), and very high critical field (> 
50 T), MgB2 has other properties that make it partic- 
ularly suitable for studying BCS superconductivity in 
some very interesting limits. One such property is the 
two gap feature that is readily observed in high quality 
crystals. Other examples include the low spin-orbit in- 
teraction and this material or the sensitivity to disorder 
that can be studied over a wider range of disorder due to 
the initial high transition temperature. 

It is the disorder effect that is the subject of the present 
paper, where samples have been fabricated in thin film 
form. These are obtained using a variant of the standard 
deposition with a post-anneal method in which the an- 
nealing is not carried through to completion, thus yield- 
ing an interesting new class of disordered (possibly amor- 
phous) superconductors with relatively high transition 
temperatures. This approach, with a controlled Mg diffu- 
sion process, yields a novel method for the fabrication of 
MgB2 /MgO multilayers. The multilayers appear to ex- 



hibit quasi-two-dimensional superconductivity and, with 
further optimization, may provide a new model system 
for studies of the effects of dimensionality on supercon- 
ductivity. 

II. SAMPLES 

As it is well known, two factors complicate the fabrica- 
tion of MgB2 films. These are the high vapor pressure of 
the magnesium at low temperature and the high sensitiv- 
ity of magnesium to oxidation, requiring very low oxygen 
partial pressures in the deposition system. In order to 
compensate for a possible loss of Mg in the subsequent 
in-situ annealing process, excess Mg was provided during 
deposition, following the approach of most groups using 
Pulsed Laser Deposition (PLD)^*^'^ This was achieved 
by successive deposition of bilayers of MgB2/Mg on a 
room temperature AI2O3 substrate. The subsequent an- 
nealing temperature was kept low in order to obtain dis- 
ordered films, with lower temperatures (of order ~ 600 
°C) yielding the most disordered films with a finite su- 
perconducting temperature. Subsequent cooldown in the 
presence of a finite oxygen partial pressure yields a novel 
multilayer structure which is spontaneously formed with 
MgB2/MgO being the unit bilayer. A schematic repre- 
sentation of our samples is shown in Fig. [H We found 
that influencing the properties of the superconducting 
layers in such a layered system was not difficult. 

Films reported in this study were grown by PLD in a 
stainless steel vacuum chamber with a background pres- 
sure of 10~^ Torr. A KrF excimer laser produces a 248 
nm wavelength beam with typical pulse lengths of 20-30 
ns. A rectangular mask shapes the beam, and a variable 
attenuator permits variation of the pulse energy. The 
energy density on the targets was kept at approximately 
2.3 J/cm^. Before each run the targets were pre- ablated 
for one minute at 20 Hertz each. 

The multilayers were grown on AI2O3 R-plane sub- 
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FIG. 1: A schematic representation of a typical sample. Mul- 
tilayers were always started with a layer of MgB2 and always 
terminated with a layer of Mg. The magnesium layers were 
found to be fully oxidized after being exposed to air. 



strates, which were annealed in an oxygen atmosphere 
at 1100 °C for one hour to create a high quality sur- 
face. Substrates were glued onto the stainless steel sam- 
ple holder with silver paste and cured for ten minutes at 
150 °C before the chamber was evacuated. The samples 
were grown by alternating between targets of sintered 
MgB2 and metallic Mg; each bilayer was deposited by 
ablating first MgB2 and then Mg. The last layer in 
the structure was made thinner than the others to en- 
able good electrical connection to the conducting layer. 
The properties of the samples, for example Tc, and the 
coupling between superconducting layers, could be influ- 
enced by varying the thicknesses of the MgB2 and MgO 
layers. A shadow mask was used to create a standard 
transport-Hall bridge pattern in order to facilitate elec- 
trical transport measurements. 

The MgB2/Mg layers were deposited at room temper- 
ature and subsequently annealed in- situ in 100 mTorr 
flowing argon at 600 °C for one hour. The samples were 
then cooled down and removed from the chamber, and 
the excess Mg oxidized. Care was taken to minimize 
the time between fabrication and measurement, to avoid 
degradation of the properties because of the exposure to 
air. 

Both low-angle (reflectivity) and high-angle X-Ray 
diffraction (XRD) experiments were performed on a 
Philips X'pert thin film diffractometer. For the reflec- 
tivity measurements unpatterned samples were used to 
maximize the signal. Electrical transport properties were 
measured with a Quantum Design PPMS system, in a 
magnetic field up to nine Tesla. Transport measurements 
at temperatures below 1 K were carried out in an Oxford 
400 dilution refrigerator. The linear sheet resistance was 
measured using standard four-point lockin techniques at 
frequencies from 3 Hz to 17 Hz, with excitation currents 
kept below 1 nA. Transport measurements results will be 
discussed in the next section. The samples reported here 
are listed in Table [J along with some of their properties. 

Figure [2] shows a typical low-angle x-ray diffraction 
measurement on one of the samples. This particular 
sample was deposited using 360 pulses of MgB2 and 120 
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FIG. 2: a) X-ray reflectivity measurement of sample A. Two 
periods are visible, indicative of a sample with two different 
types of interfaces. The best simulation (dotted line) was 
obtained when only 7 bilayers were included in the model. 
The densities of the layers were close to the theoretical values 
of MgB2 and MgO and the roughness was increased linearly 
from the substrate to the top surface. The simulation is offset 
for clarity, b) High angle XRD measurement of sample A. 
The main peak is due to the AI2O3 substrate. Arrows point 
to where MgB2 lines are expected. 



pulses of Mg for every bilayer (sample A in Table [I]) , re- 
peating this sequence eight times. Apparent in Fig. [2] are 
oscillations of two distinct periods. We interpret these 
oscillations in terms of a superstructure consisting of re- 
peating bilayers. The oscillations with the large period 
come from the layered pairs, and those with the short 
period come from the entire film. 

From the period of the large peaks, we estimate the 
thickness of a bilayer to be approximately 11 nm. Since 
the film was created by alternately pulsing between 
MgB2 and Mg eight times, one would expect that there 
are eight repeats of the bilayers. However, the best simu- 
lation (the dotted line in Fig. [2^) was obtained using only 
seven such layers, resulting in a total film thickness of 77 
nm. We believe this discrepancy in the total number of 
layers is due to a roughening of the superstructure that 
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TABLE I: Characteristics of samples reported in this paper. 
Multilayers consist of 8 bilayers of MgB2 /MgO. 



Sample 


MgB2 


Mg 


d-MgB2 


d-Mg 


d-film 


Rn(40K) 


Tc 




pulses 


pulses 


[nm] 


[nm] 


[nm] 


M 


[K] 


A 


360 


120 


8.0 


3.0 


77 


91.3 


24 


B 


360 


600 


8.0 


15 


185 


46 


17 


C 


240 


600 


5.3 


15 


163 


786 


9 


D 


220 


600 


4.9 


15 


160 


1363 


7 


E 


180 


600 


4.0 


15 


153 


420000 





occurs with increasing thickness, causing the top layer 
to be destroyed and not show up in reflectivity measure- 
ments. SEM measurements support this roughening con- 
jecture and show a relatively flat template on which large 
(presumably boron^) particles are present. The simula- 
tion of the reflectivity measurements also tells us that the 
densities of the materials are very close to the theoretical 
values expected for MgB2 and MgO, indicating that the 
magnesium in the MgO layers has fully oxidized. To get 
a good simulation for the higher angles (> 1°) the RMS 
roughness was increased linearly from the substrate to 
the top surface. 

High angle XRD measurements show the substrate 
peak and some very broad low intensity peaks that can 
be attributed to MgO, as is shown in Fig. [2)3. The lack 
of strong MgB2 peaks confirms that the MgB2 is highly 
disordered and possibly even amorphous. 



III. NORMAL STATE TRANSPORT AND THE 
SUPERCONDUCTING TRANSITION 

The normal state resistance is a first good indicator 
of the disorder in our samples. For high quality single 
crystals the resistivity has been found to drop with tem- 
perature, exhibiting a temperature dependence Ap oc 
with a > 2.5 and resistivity ratio exceeding 257*i^^^ 
Similar behavior was also found for high quality, dense 
wires of MgB2.^ Quality thin films, usually grown epitax- 
ially, exhibit similar normal state properties.^ However, 
upon introduction of disorder, single crystals show a di- 
minished resistivity ratio. For carbon-doped crystals this 
ratio was reduced to order unity, while the superconduct- 
ing transition temperature decreased monotonically with 
increasing carbon content 

Inspecting the normal state resistance of our films in 
Fig. [31 we observe that the samples all show very low re- 
sistivity ratios, similar to disordered bulk materials. Note 
that we have calculated the two-dimensional resistivity, 
i?n, for the full thickness of the sample and not per layer. 
The resistance per square for films A (not shown) and B 
is low even though the resistivity is of order 10~^ l]-cm, 
which no doubt explains the almost flat temperature de- 
pendence of the resistance. The transition temperature 
Tc of these two samples is reduced by a factor of 2-3 from 



that of bulk, high-quality crystals Tc. Samples C and D 
are more anisotropic and consequently show larger resis- 
tance. These two samples also show a resistance mini- 
mum at ~150 K, signaling the onset of localization pro- 
cesses. Indeed, a slight change in the thickness of the 
MgB2 superconducting layers, from 4.9 to 4.0 nm, drives 
the sample from being metallic and then superconduct- 
ing (sample D) to being insulating (sample E), despite 
the same MgO spacer layers and same total number of 
layers. 
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FIG. 3: Normal state resistivity of samples B, C, D (top) 
and E (bottom). Note the different scale for the two panels. 
For the top panel we also show the transitions in a magnetic 
field. These are discussed in the text. 

In Fig. [3^1 we plotted the normal state resistance in 
magnetic fields between zero and 9 T, visible on the left 
side of the graph. While we examine the superconducting 
transitions in a magnetic field in greater detail in the next 
section, we note here the rather simple behavior of the 
transitions in low fields as being simply superconducting 
pair-amplitude controlled. This is clearly seen by the 
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fact that the normal state resistance seems to continue 
smoothly to lower temperatures with increasing field. 

Figure 2] concentrates on the zero-field resistive tran- 
sitions of the representative samples listed in Table [D 
Fig. [4^ shows the transitions for various thicknesses of the 
MgB2 layers, with a fixed MgO thickness, while Fig. [4)3 
shows the transitions at various MgO thicknesses, with a 
fixed MgB2 thickness. Two prominent features of these 
data are that the normal state sheet resistance is roughly 
independent of temperature and that is reduced from 
the bulk value of 39 K, both as expected for disordered 
superconductors. Nevertheless, in absolute terms the Tc's 
are distinctly high for a disordered conventional super- 
conductor. Also evident from Fig. [4^ is that the normal 
state sheet resistance increases and Tc decreases as the 
MgB2 layer is made thinner. This may be attributed to 
the stronger two-dimensionality and thus weaker super- 
conductivity of the MgB2 layers, unable to strengthen 
due to the thick, intervening MgO layers. This interpre- 
tation is reinforced by Fig.[3]3, showing that Tc decreases 
as the thickness of the MgO layers increases, further de- 
coupling the MgB2 layers. 



IV. UPPER CRITICAL FIELD 

The best way to study the coupling between layers 
is by using a magnetic field. Indeed, magnetic prop- 
erties of high-temperature superconductors, especially 
in the vortex state, revealed much of their anisotropic 
character. ^^^^^ This reality led White et al. to propose 
that artificial multilayers of MoGe/Ge could be used as 
a tunable model system for the study of the magnetic 
properties of layered superconductors 

The anisotropy of the multilayers is demonstrated 
through the angular dependence of the resistance near 
and below Tc. Figure [5] shows the strong anisotropy of 
sample A, as the temperature is lowered below the zero- 
field Tc at an applied field of 9 T. This sample, with the 
highest Tc in Table [H has relatively strong anisotropy 
demonstrated by the fact that a few degrees below Tc 
the parallel field resistance is zero while the perpendicu- 
lar one is a significant fraction of the normal state resis- 
tance. 

Typical results for the magnetic field dependence of the 
resistive transition and the associated upper critical field 
curves Hc2{T) (perpendicular and parallel to the layers) 
are shown in Fig. [6] for sample B. A set of resistive transi- 
tions in perpendicular fields up to 9 T shows broad transi- 
tions that are roughly parallel with increasing field. From 
this data we can determine the in-plane coherence length. 
Choosing a mid-point of the transition, the perpendicular 
critical field (see inset in Fig. [6^) shows a linear behav- 
ior with temperature, near Tc. From the slope of the 
critical field near Tc we calculate the zero-temperature 
coherence length using: -Tc{dHc2±/dT) = <l>o/(27r^^). 
We find that the in-plane zero-temperature Ginzburg- 
Landau (GL) coherence length ^ 50A, about 50% 
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FIG. 4: a) The zero-field resistance for different MgB 2 su- 
perconducting layer thicknesses, keeping the MgO separating 
layer constant at 600 pulses. Note the logarithmic scale. From 
bottom to top the curves of samples B,C,D and E are plot- 
ted. The number of pulses for the MgB2 layer is given for each 
curve relative to the 360 used for sample B. b) Comparison 
of samples A and B with 360 pulses of MgB2 but different 
number of Mg pulses, resulting in a different spacing between 
the superconducting layers. The more decoupled sample (A) 
shows the greater reduction in Tc, down to 17 K, compared 
to 24 K for the other sample. 



lower than the in-plane coherence length found for high 
quality bulk samples of MgB2.^ Resistive transitions in a 
parallel field also shows parallel curves, with Tc decreas- 
ing as the field increases. Again we choose a mid-point 
of the transition and plot the parallel critical field in the 
inset of Fig. [6)3. 

For this parallel-field case we note that two anisotropy 
effects need to be considered. The first is the inherent 
anisotropy of the MgB2 material which is of order 6 for 
high quality crystals but was found to be substantially 
reduced when disorder is introduced. Indeed for dis- 
ordered crystalsi^ or films ^ the inherent anisotropy be- 
tween the ah and c directions of the MgB2 can be reduced 
to unity. The second effect is the anisotropy due to the 
layered nature of our sample. The very steep increase 
of the parallel field near Tc prevents us from determining 
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FIG. 5: Angular dependence of the resistance at H =9 T for 
sample A. The strong anisotropy is evident in the 18 K data. 



the coherence length in that direction. We conclude that 
this coherence length has to be much smaller than the 
single MgB2 layer thickness, and the mass-anisotropy ra- 
tio cx {iin/iout)'^ >> 1. The large anisotropy observed is 
therefore a direct result of the layered nature of our sam- 
ple and not the inherent anisotropy of the MgB2 layers. 

We next turn to sample D which is more disordered 
than sample B based on its normal state resistivity. Re- 
sistive transitions in a field for sample D are shown in Fig- 
ure [71 Compared with the transitions for sample B, those 
of sample D are significantly broader. However, since the 
transitions are still roughly parallel with increasing field, 
we can again determine in-plane coherence length from 
the slope of Hc2{T) versus Tc. Choosing again a mid- 
point of the transition, the perpendicular critical field 
shows a linear behavior near Tc, and from the slope near 
Tc we calculate ~ 60A. 

For the parallel field transitions, we also choose a mid- 
point of the transition for the calculation of coherence 
lengths. The inset in Fig. [7)3 shows a strong curvature 
behavior below Tc characteristic of two-dimensional sys- 
tems. The behavior appears less steep than for sample B 
and thus we attempt to analyze the parallel critical field 
behavior as a crossover from three to two dimensions, for 
which the critical field near Tc has the temperature de- 
pendence Hc2\\ oc (Tc-T)^/^. In the limit as T ^ Tc, one 
expects the three-dimensional result Hc2 oc (Tc — T) . For 
a multilayered system such behavior is a signature of de- 
coupled layers. Indeed, close inspection of the behavior 
near Tc suggests that below ~ 2 T a straight line can be 
fit to the data with a slope of ~ 2.4 T/K yielding a coher- 
ence length perpendicular to the planes of (^out ~ 13A, 
smaller than the superconducting layer thickness. The 
ratio of the two coherence lengths is ~ 4.6 (implying 
a mass-anisotropy value of ~ 20) which, as explained 




FIG. 6: a) Resistive transitions of sample B in a perpen- 
dicular magnetic field. The fields range from to 9 T, in 
increments of 1 T. The inset shows the dependence on tran- 
sition temperature of the in-plane critical magnetic field, Hc2 
(taken at 50% of the transition). The coherence length was 
calculated to be 51 A. b) Resistive transitions in a paral- 
lel magnetic field. The inset shows the critical magnetic field 
plotted against the transition temperature at 50% of the tran- 
sition. 



above, should be viewed as the anisotropy of the mul- 
tilayered system rather than the intrinsic anisotropy of 
MgB2. 

The pure two-dimensional regime of the parallel crit- 
ical field can be obtained by plotting the square of the 
parallel critical field versus temperature. For the data 
in Fig. [7)3, a linear regime is obtained below ~ 14 K. 
In this regime we can use the expression^'^ Hc2\\{0) = 

^o/{27T^abd)^ where d is some effective thickness of order 
the superconducting layer thickness. Using the data, we 
obtain d ~ 33A, a value that is close to the thickness of 
the superconducting layer (49 A, see table [I]). We note 
however that the data do not extrapolate to the mea- 
sured Tc. One possible reason for the discrepancy is that 
since the two-dimensional regime occurs when the inter- 
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FIG. 7: a) Resistive transitions of sample D in a perpendic- 
ular magnetic field. As in Figure 6, the fields range from 
to 9 T, in increments of 1 T. The inset shows the dependence 
on transition temperature of the in-plane critical magnetic 
field, Hc2 (taken at 50% of the transition). The coherence 
length was calculated to be 58 A. b) Resistive transitions in a 
parallel magnetic field. The inset shows the critical magnetic 
field plotted against the transition temperature at 50% of the 
transition. 

layer coherence length decreases below the bilayer thick- 
ness the lower Tc may point to the interaction-reduced 
Tc typical for thin superconducting films 

V. LOW TEMPERATURES AND HIGH FIELDS 

We turn now to the analysis of behavior of the multi- 
layers at very low temperatures and high magnetic fields. 
Figure [8] shows a series of isotherms of the resistance as a 
function of magnetic field, for sample D. The isotherms 
each increase in resistance above some characteristic field, 
and appear to saturate at high field. The value of this 
saturation resistance increases with decreasing tempera- 
ture, which is an indication of increased localization ef- 
fects with lower temperature. Another interesting fea- 
ture of the data is that the isotherms all cross at a point 
(Hc^Rc)- In general, such a crossing is a signature of 



a zero-temperature quantum phase transition for which 
the magnetic field is a tuning parameter. 

Indeed, disordered superconducting films are expected 
to undergo a phase transition to an insulating state in 
the limit of zero temperature, tuned by either a magnetic 
field or by strong disorder. In the case of the field-tuned 
transition, one finds that the linear resistance decreases 
with decreasing temperature when the field is below a 
critical field i^c, and increases with decreasing tempera- 
ture when the field is above the critical field. However, 
the insulating behavior observed in our films above the 
crossing point is very weak, and in fact similar to the 
behavior observed in high temperature superconductors 
as was previously described by Steiner et al^ 
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FIG. 8: Resistance of sample D as a function of perpendicular 
magnetic field for 400 mK, 200 mK, 120 mK and 60 mK. The 
curves show a crossing point at Re — 10.42 T and Rc — 
2350 n/D. 

The crossing point observed in Fig. [8] also relates to the 
determination of the upper critical field for this sample. 
Note that in the previous section, using the mid-point 
of the transition we found an almost linear dependence 
of Hc2{T) all the way to zero temperature, extrapolat- 
ing to Hc2{0) ^ 9.2T. This value is significantly less 
than the crossing point field of He = 10.42 T. If, in- 
stead of using the mid-point we use the actual zero re- 
sistance point of the sample to determine iJc2, we get 
an initial slope of 0.68 T/K (implying an in-plane co- 
herence length of 87 A), which appears to become very 
steep below ~ 1.5 K. Similar behavior has been observed 
previously for the resistance of high-Tc superconductors 
with large anisotropy.^^ This behavior is consistent with 
the high field crossing point and, since MgB2 is believed 
to be a conventional BCS superconductor, this suggests 
that the anomalous behavior of Hc2 is a property of 
the quasi-two-dimensional layered nature of the mate- 
rials rather than the mechanism for high-Tc. Indeed 
Kotliar and Varma^ have suggested that the anoma- 
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lous behavior of Hc2 near zero temperature observed in 
some high-temperature superconductors is a consequence 
of the proximity to a zero-temperature quantum critical 
point, which is the zero-temperature critical point at the 
second-order end point of the first-order melting line of 
the vortex lattice. Shimshoni et al?^ made a similar con- 
nection to the superconductor-insulator transition in thin 
disordered films of MoGe^— 



VI. CONCLUSIONS 

In this paper we demonstrated a new method for pro- 
ducing multilayers of MgB2 with MgO as the insulating 
spacing layer. We have shown that the anisotropy of the 
multilayers can be controlled with the spacing of the insu- 
lating layer, while Tc is a function of the disorder in the 
individual layers and the coupling between layers. We 



have also shown that these materials exhibit all of the 
characteristic features of disordered superconductors, in- 
cluding a signature of a superconductor-insulator transi- 
tion. We hope that these multilayer structures could be 
used as a model system for a future study of disordered 
quasi- two-dimensional superconductors. 
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